Quantum dynamics of cavity assisted photo-association of Bose-Einstein condensed 

atoms 
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We explore the quantum dynamics of photoassociation of Bose-Einstein condensed atoms into 
molecules using an optical cavity field. Inside of an optical resonator, photoassociation of quantum 
degenerate atoms involves the interaction of three coupled quantum fields for the atoms, molecules, 
and the photons. The feedback created by a high-Q optical cavity causes the cavity field to become 
a dynamical quantity whose behavior is linked in a nonlinear manner to the atoms inside and where 
vacuum fluctuations have a more important role than in free space. We develop and compare several 
methods for calculating the dynamics of the atom-molecule conversion process with a coherently 
driven cavity field. We first introduce an alternate operator representation for the Hamiltonian 
from which we derive an improved form of mean field theory and an approximate solution of the 
Heisenberg-Langevin (HL) equations that properly accounts for quantum noise in the cavity field. 
It is shown that our improved mean field theory corrects several deficiencies in traditional mean 
field theory based on expectation values of annihilation/creation operators. Also, we show by direct 
comparison to numerical solutions of the density matrix equations that our approximate quantum 
solution of HL equations gives an accurate description of weakly or undriven cavities where mean 
field theories break down. 



I. INTRODUCTION 

In the last decade, there has been considerable interest 
in creating ultra-cold quantum degenerate molecular sys- 
tems because of the potential for improved understand- 
ing of molecular physics and interactions, exploring new 
types of many-body systems such as condensates of dipo- 
lar molecules and the BCS-BEC crossover, and the gen- 
eration of entangled atoms by controlled dissociation of 
molecules Two techniques, magnetically tunable Fes- 
hbach resonances and photoassociation, have been devel- 
oped to create ultra-cold molecules that start first from 
laser cooled atoms and then induce controlled chemical 
bonding between the atoms. Feshbach resonances are 
the most widely used and have been successfully applied 
by numerous research groups to create molecular dimers 
starting from either a Bose-Einstein condensate (BEC) 
@, 1 II 1 or a Fermi gas @, 0, B H- This work cul- 
minated in the formation of a molecular Bose-Einstein 
condensate (MBEC) HSHIl. 

Besides Feshbach resonances, experiments have 
demonstrated that two-photon Raman photoassocia- 
tion can also be used to create ultra-cold molecules 
[l2l [H, [3, EH- Two-photon photoassociation has the 
added benefit that the frequency difference between 
the two optical fields can be used to select a particu- 
lar rotational-vibrational state including the rotational- 
vibrational ground state [l|| [lj]]- This gives photoassoci- 
ation an advantage over Feshbach resonances since Fesh- 
bach molecules are often very weakly bound in high en- 
ergy vibrational states that quickly decay to lower lying 
vibrational states via inelastic collisions 

HE [Hi Re- 
cent work has has used two-photon Raman transitions to 
create rovibrational ground state molecules from weakly 
bound Feshbach molecules ISM SjJ ■ 



Here we address the problem of two-photon Raman 
photoassociation of an atomic BEC inside of an optical 
cavity that is coherently driven. Due to the cavity, pho- 
tons circulate and interact with the atoms and molecules 
many times before finally exiting in a manner analogous 
to a feedback loop. This feedback amplifies the back 
action of the atoms and molecules onto the cavity field 
causing the light to now become a dynamical part of the 
process. In our particular case, we assume that one of 
the optical fields used to induce the atom-molecule con- 
version is a quantized mode of a driven Fabrey-Perot res- 
onator while the other field does not correspond to a cav- 
ity mode but is rather a laser directed transverse to the 
cavity with sufficient intensity to be treated as a 'classi- 
cal' undepleted pump. Our model therefore involves the 
interaction of four particles: two atoms are 'destroyed' 
and a molecule and cavity photon are 'created' and vice 
versa. Consequently, the atom-molecule-cavity photon 
interaction is analogous to susceptibility in nonlin- 
ear optics. This is different from molecule formation via 
a Feshbach resonance or free space photoassociation with 
undepleted classical lasers where the conversion only in- 
volves two quantum fields: atoms and molecules and is 
the matter- wave analog of second harmonic generation of 
photons with a x^ 2 ) susceptibility. Coherent photoasso- 
ciation inside of a cavity therefore offers the prospect of 
novel nonlinear dynamics between the atomic, molecular, 
and cavity fields and the possibility of enhanced control 
over the atom-molecule conversion process. 

In an earlier work we analyzed the mean field dynamics 
and steady state behavior of cavity assisted photoassoci- 
ation 22] while here we extend that work to study the 
role of quantum fluctuations on the dynamics. We an- 
alyze the quantum dynamics for the atomic, molecular, 
and cavity fields by several methods. First, an alter- 
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nate operator representation for the Hamiltonian is intro- 
duced that has an algebra analogous to angular momen- 
tum. These new operators allow us to derive from the 
Heisenberg-Langevin equations both an improved form 
of mean field theory and an approximate solution for 
the quantum dynamics that treats the atom-molecule 
populations classically but the cavity field and atom- 
molecule-photon coherences fully quantum mechanically. 
The improved mean field theory incorporates quantum 
correlations between the atom and molecules. It also in- 
cludes a contribution due to vacuum fluctuations of the 
atomic and molecular fields that corrects a deficiency 
in traditional mean field theory, which fails to predict 
atom-molecule Rabi oscillations for resonant transitions 
because the solution approaches and becomes stuck at 
an unstable equilibrium point. We also compare the im- 
proved mean field theory with our approximate quan- 
tum solution and direct integration of the density matrix 
equations in the case of small photon and atom numbers. 
It is shown that our approximate quantum solution pro- 
vides a more accurate description in the case that the 
cavity driving is weak or absent in comparison to mean 
field theories since in this case the dynamics are initiated 
by the vacuum fluctuations of the cavity field. 

Before proceeding we note that only a few papers have 
previously considered photoassociation inside of a cavity 
[23L f2^| . However, unlike our model, theirs was based on 
single photon photoassociation, which is impractical for 
observing coherent atom-molecule dynamics because the 
molecules created are in electronic excited states and can 
rapidly decay due to spontaneous emission. The authors 
of Ref. [H 111] employed the positive-P distribution [H 
to analyze the quantum dynamics of the three coupled 
fields. The positive-P distribution has a tendency to be- 
come numerically unstable for long times in highly non- 
linear quantum optical systems as the stochastic trajec- 
tories begin to sample unphysical regions of phase space 
and must be stabilized by proper choice of a stochastic 
gauge, which is often difficult to properly determine [26| . 
In fact, it has been shown that the equations presented 
in the earlier work [23|, [24| are numerically stable only 
for a limited range of parameters and that the photon 
blockade effect predicted in that work does not in gen- 
eral occur [271]. 

The paper is organized as follows: In section II, we 
present our model for cavity assisted photoassociation 
and the density matrix equations of motion. Additional 
details for the physical model can be found in Ref. [22| . 
In section III, we derive our improved mean field theory 
and an approximate solution of the Heisenberg-Langevin 
equations for the dynamics based on a pseudo-angular 
momentum operator representation. In section IV, we 
compare numerical results for the density matrix, approx- 
imate Heisenberg-Langevin equations, and mean field 
theories. 



II. MODEL 

Our starting point is a BEC of atoms inside of an op- 
tical cavity, as depicted in Fig. [TJ The atoms as well as 
the molecules formed from them can be trapped inside of 
the cavity using a far-off resonant optical trap similar to 
what has been recently demonstrated with single atoms 
in a cavity 28]. At temperatures T m 0, we can assume 
that all of the atoms are in the ground state of the trap- 
ping potential with wave function ij) a (r). Additionally, 
the atoms are assumed to have all been prepared in the 
same hyperfine state denoted by \a). Pairs of atoms in \a) 
are coupled to electronically excited molecular states \I V ) , 
where v denotes the vibrational state of the molecule, via 
a pump laser with Rabi frequency fli and frequency uii. 
The pump is treated as a large amplitude undepleted 
source and therefore changes in fi; due to absorption or 
stimulated emission are neglected. 




FIG. 1: (Color online) Schematic diagram showing the system 
under consideration. The kets \a), \b), and \I V ) denote pairs of 
atoms, electronic ground state molecules, and electronically 
excited molecules, respectively, e is the rate at which the 
cavity is coherently driven by an external laser while 7 is the 
decay rate for photons in the cavity, fli is the pump laser that 
drives the \a) — > \I V ) transition. 

The excited molecular states are coupled to molecules 
in their electronic and vibrational ground state, \b), via 
a single cavity mode. The ground state wave function 
for the center of mass of the molecules is denoted by 
^b(x). Emission of a photon into the cavity mode takes 
a molecule from an excited state to its electronic ground 
state. Coupling to a single mode can be achieved by 
insuring that only a single cavity mode is close to two- 
photon resonance for the atom-molecule Raman transi- 
tion and by positioning the atoms and molecules around 
an antinode of the cavity field. The discrete mode struc- 
ture of the cavity allows one to select a particular vi- 
brational state in the electronic ground-state manifold of 
the molecules provided the cavity line width, 7, is less 
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than the vibrational level spacing, which can be as high 
as 1GHz [l6|]. This would imply a cavity Q-factor of 
Q 10 6 , which has already been achieved with individ- 
ual atoms trapped inside of a Fabrey-Perot resonator (28). 
The cavity field frequency is oj c and the vacuum Rabi fre- 
quency for the |J„) — * | b) transition is g ca v 

The internal energies of states and \b) relative to 
pairs of atoms in \a) are uj v and Auj < 0, respectively. We 
assume that the detuning between the excited states and 
the pump and cavity fields satisfy, A v = u>„ — u>i ~ — 
Aw)-w c 3> |g C a-u|, 1, 7i/ where 7" 1 is the lifetime of \l v ) 
due to spontaneous emission. Under these conditions the 
excited states can be adiabatically eliminated, leading to 
two-photon Raman transitions between \a) and \b) with 
the resulting Hamiltonian for the atom-cavity-molecule 
system, 



H, 



hsPb 



+ i — (a t2 6e - d 2 P ^ ) + h\ a a 



t 2 a 2 - 



hxbb^b 2 . 
(1) 

Here a, 6, and e are bosonic annihilation operators for 
atoms, ground state molecules, and cavity photons, re- 
spectively. Moreover, the molecular and photon op- 
erators have been written in a rotating frame, b — > 
bexp[+i(u) c — u>t)t] and e — > eexp[— iu) c t\, to remove all 
time dependence from the interaction term. The two- 
photon detuning is defined as 5 — (Alu — (cjj — u> c )). 

The terms proportional to Xa and Xb represent the 
two body interactions between pairs of atoms and pairs 
of molecules, respectively. Interactions involving an atom 
interacting with a molecule can be written as XabN a N{, = 

jXab (n 2 -N 2 - 47V 2 ) where N = N a + 2N b is the total 

number operator and N a = cva and N b = tfb. Since the 
total number of particles is conserved, [H am , N] — 0, Xab 
can be absorbed into a redefinition of Xa and Xb- 

The dynamics of the cavity are described by two com- 
peting processes. The first process is cavity decay, which 
can be treated using the standard Born-Markov master 
equation for the density operator [25| . 



dp 



\damping 



j- (2e / oe t - e f ep - /oe f e) (2) 



In addition to this, the cavity is coherently driven by an 
external laser described by the following Hamiltonian, 



(3) 



where we have assumed that the driving laser is resonant 
with the cavity mode. The complete equation of motion 
for the density operator is then given by, 



H am , p] + 7: (2epe f - e f ep - pe 



dt~ih lpump 

(4) 

The two-photon Rabi frequency is g 
-i J d z xi)l(x)[g cav u{x)W l {x)Y JU Il, u IbMl 'K]^a{x). I l>v 
are the Frank-Condon factors for the \l — a,b) — > |/„) 
transitions [2i| and since typically |f2;|/|A„| and 



Ia,u -C 1, 5/7 <C 1 even if the cavity is in the strong 
coupling regime \g cav \ > 7 

We represent the density operator in the basis of eigen- 
states of N a , and N e = e^e, |n , n e ). Because n a + 2rib = 
N is a constant of motion, the molecule number, tij,, is 
completely determined by n a . The density matrix in this 
basis, (n' a , n' e \p\n ai n e ), is then unwrapped into a column 
vector, /?, so that Eq. 3] can be written as a matrix equa- 
tion, dp/ dt = Mp, which can now be integrated using a 
first order Euler method 1301] . 



p(t) = lim (7 + (t/k)M) k p{0) 



(5) 



This method has the advantage of using less memory 
than higher order ODE solvers but is still limited to small 
numbers of atoms and photons. In the next section we de- 
velop approximate equations of motion that incorporate 
quantum effects but can deal with much larger experi- 
mentally realistic numbers of atoms (~ 10 3 — 10 6 ) and 
photons. 



III. PSEUDO-ANGULAR MOMENTUM 
DESCRIPTION 

Here we develop a representation of the model devel- 
oped in the previous section using pseudo- angular mo- 
mentum operators that simplify the form of the Hamil- 
tonian and use this representation to derive Heisenberg- 
Langevin equations for the system. We then obtain mean 
field equations and approximate equations for the quan- 
tum dynamics from the Heisenberg-Langevin equations. 

For an initial state that is an eigenstate of N = N a + 
2Nb, the solution of Eq. 2] will at all later times continue 
to be an eigenstate of N with the same eigenvalue N. In 
this case, we can introduce new operators [3l| 



L, 



V2 (tf&b + tfdaj 



7V3/2 



(6) 



s/2[tftfb-tft 
<L + -L_)= V — '- (7) 



iN 3 / 2 



L z = 



- tic 
N 



(8) 



and V + = These operators have the following com- 
mutation relations, 

L±,£„] = ±jjl± (9) 
r~ . 1 3 /- \ f. \ 2 



±(L.-i) (4 + 1/3) -A (10) 



The commutation relations 



L-t-, L, 



are of the same form 



_ as for angular momentum but the equivalence to angular 



We 



momentum is ruined by the commutator 

therefore refer to them as pseudo-angular momentum op 
erators. 
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Equation [T] can be rewritten in terms of these new op- 
erators as 



H am = ^N5L z + ^N 2 X L 2 z + ihg (E 



3/2 



L+e-L-eJ 



(11) 

where X = Xa + W 4 and 8 = S + N{-2 Xa {l - l/N) + 
( Xfe /2)(l-2/AT)). 

It is well known that Born-Markov density ma- 
trix equations such as Eq. U are fully equivalent to 
Heisenberg-Langevin equations with Markov noise oper- 
ators in the Heisenberg picture (25[, which in this case 
have the form 



dt 



d L 

Jt* 



iN X (L z L_+L_L; 



■9(2N) 



N 



1/2 



AJ V 2 

-g{2N) 1 ' 2 (L + e + L 



L z -1) (L, + l/3je (12) 
n (13) 



_e 



dt 



2 6 3 V 2 



3/2 



L_+F(t) 



(14) 



Here F(f) is a Markov noise operator for the fluctua- 
tions of the electromagnetic reservoir coupled to the cav- 
ity mode via the mirrors with zero mean (F(t)) = and 
the two-time correlations (F(t)F(f)) = (Ft(t)F(i' )) = 
and (F(f)Ft(f)> = -y6(t-t'). 

Because the Heisenberg-Langevin are nonlinear op- 
erator equations, they cannot be solved exactly. The 
simplest approximation is that of mean field theory, 
which replaces all operators with their c-number ex- 
pectation values and factorizes products of operators, 
L + e — ► (L +)(£.} thereby ignoring higher order correla- 
tions and non-commutativity of operators. This yields 
three nonlinear c-number differential equations, 



dt 



dt 



iSL_ + i2N X L z L_ + g(2N)- 1/2 e 
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N 



1/2 



dt 



4/ V 2 

-g(2Nfl 2 (L +e + L_e*) 

7 (N\ 3/ \ 



(L,-l)(£, + l/3)e(15) 
(16) 
(17) 



where the lack of a hat denotes a c-number expectation 
value, A = (A). These equations are unaffected by quan- 
tum noise since (F(t)) = and according to mean field 
theory, (ee^) = (e^e) = e*e. 

We note that this form of the mean field approximation 
is different from the traditional manifestation with Bose- 
Einstein condensates or light where annihilation/creation 
operators for the fields are replaced with c-number expec- 
tation values. This form was used in our previous work 
[2^ | and yielded the equations 



a = —2i Xa a* a 2 + ga*be 



(18) 



b = -i5b~2i Xb b 2 b* - -a 2 e* 
2 

7 9 2,* 

e = e e a 

2 2 



(19) 
(20) 



where a = (a), b — (b), and again e = (e). In the case 
that e(t) is treated as a time independent constant, Eqs. 
[TBI and are equivalent to earlier mean field studies 
of atom-molecule dynamics with Feshbach resonances or 
photoassociation [l], 0, [32|, [33|. By contrast Eqs. 
I15H171 which we refer to as pseudo-angular momentum 
mean field theory (PAMMF), represent a higher order 
approximation for the atomic and molecular fields than 
Eqs. [T8"ll2T)l which we refer to as amplitude mean field 
(AMF) theory. The PAMMF equations automatically 
incorporate lowest order quantum correlations between 
the atomic and molecular field operators as a conse- 
quence of the pseudo-angular momentum representation, 
L z = (2ttb-tfa)/N and L_ = V2(tf da) /N 3 ? 2 . Addi- 
tionally, Eq. [TS] includes the vacuum fluctuations of the 
matter fields resulting from the last term in the com- 
mutator In the next section we consider the 
different predictions for the dynamics made by PAMMF 
and AMF. 

Mean field theory typically works well for large ampli- 
tude quantum fields. However, initially the cavity field 
is in the vacuum state and therefore it cannot be used 
to describe the initial short time behavior ~ 7 ■ Addi- 
tionally, when the cavity is not driven e = 0, mean field 
theory completely fails to describe the dynamics, which 
are initiated by vacuum fluctuations of the cavity field. 
Therefore we develop a solution to Eqs. [T2"l[l"4"l that prop- 
erly includes photon noise. 

We note that equations [T2l and fT4l can be linearized and 
solved exactly if one replaces L z with a c-number, L z . 
In the case that the c-number is time independent, this 
approach would be equivalent to the undepleted pump 
approximation in nonlinear quantum optics and the Bo- 
goliubov method for weakly interacting BEC's. We adopt 
an approach where the c-number L z is instead treated as 
a dynamical variable whose value is given by the expec- 
tation value of Eq. 1131 First we express the equations 
for and e as 



—X = MX + S 
dt 



where X T = ( L_ , e ) and S T = ( 0, e + F(t) ) while 



(21) 



_ , iAQ(t) -g/3(t) \ 
M ^ - I -g{N/2f/ 2 - 7 /2 



(22) 



and A£l = 6 
(3/4)(JV/2)V a (L,(t) a . 
The correlation matrix, 

C(t) = (X(t)X(tf) 



2N x L z (t) 
2L z {t)/3 - 



(L-L+) 
(eL + ) 



(23) 
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gives us information about the number of cavity photons 
and the coherence between the cavity field and the atom- 
molecule fields, which can be used to calculate L z (t). 
C(t) has the equation of motion obtained directly from 
Eq. M 



dt 



C = MC + CM* 







s*(L_) 
7 + e(e)* +e* 



(24) 

Equation [2U forms a closed set along with the equations 
for the expectation value of Eq. [5T] and the expectation 
value of L z (t) , 



dt 



L z = 



-g{2N) 1 ' 2 (<eL 4 



(L-e 



at- 



(25) 



We refer to this approach that incorporates the quan- 
tum noise of the cavity field and atom-molecule-photon 
coherences while treating the atom-molecule popula- 
tions classically as the quantum self-consistent popula- 
tion (QSCP) method. In general this QSCP method is 
valid for large N ^> 1 in the same manner as the Bogoli- 
ubov theory of weakly interacting condensates [l[ • From 
the c-number substitution L z — ► L z (t), one can see then 
that the commutation relations Eqs. [9] and [10] are only 
preserved in the limit N — > oo. In comparison to AMF 
and PAMMF equations, Eqs. [24] and [25] properly ac- 
count for quantum noise of the cavity field. The quan- 
tum noise of the photons arise from two sources in the 
equations for C(t): (ee*) — {e'e) + 1 representing the 
vacuum fluctuations of the cavity mode and the reser- 
voir noise (F(t)F^(t')) = jS(t — t'), which appears as 7 
in the equation for (ee'). 



state Nb = N/2 molecules and N e w photons. For 
resonant atom-molecule conversion (S — \a — Xb = 0), 
PAMMF and QSCP agree very closely for the first few 
Rabi oscillations while for longer times the QSCP so- 
lutions oscillate at a slightly higher frequency. For off- 
resonant transitions (6 7^ or \a, Xb 7^ 0), the agree- 
ment becomes increasingly better and they are indistin- 
guishable for \S\ 3> gVN- This indicates that quantum 
correlations between the cavity field and atom-molecule 
medium given by (eL + ) play a non-negligible role in mod- 
ifying the dynamics when full molecule conversion can 
occur. 



x 10 




IV. NUMERICAL RESULTS 

In this section we analyze numerical solutions of the 
AMF, PAMMF, QSCP, and exact density matrix equa- 
tions. In all simulations we use the initial conditions 
N a (0) = N atoms, Nb(0) = molecules, and no photons, 
N e (0) - 0. 



A. Semiclassical Limit, N ^> 1 and e 3> 7 

When the cavity mode is strongly driven, e 3> 7, the 
photon noise is negligible and atom-molecule transitions 
are dominated by stimulated absorption and emission in- 
stead spontaneous emission into the cavity mode imply- 
ing that one should be able to treat the cavity mode 
classically. Here we compare traditional mean field the- 
ory given by AMF with the improved PAMMF mean field 
theory along with the the QSCP method that incorpo- 
rates cavity noise. 

Figure 2 shows that the PAMMF mean field dynamics 
(Eqs. fT5HTT|) reproduce qualitatively the behavior of the 
the QSCP method (Eqs. [M] and |2"5"]) . which consist of 
damped Rabi oscillations that approach the final steady 



FIG. 2: (Color Online) (a) Molecule number and (b) intra- 
cavity photon number calculated using the PAMMF mean 
field equations (solid lines) and QSCP equations (dash dot 
lines) as a function of time t (in units of 7 )• In all cases 
N = 10 6 , 3/7 = 5 x 1CT 5 , and e/7 = 100. Here red lines are 
for S — Xa = Xb = while green lines are for 5/7 = 1, and 
Xa = Xb — 0. Note that in the later case (green lines), the 
two solutions are indistinguishable. 



By contrast, the AMF solutions (Eqs. [T8][20]) exhibit 
qualitatively very different dynamics from PAMMF and 
QSCP method in the resonant case as seen in Fig. 3. 
In the case 6 — Xa — Xb — 0, AMF equations dis- 
play no Rabi oscillations and instead the molecule num- 
ber grows monotonically until it reaches the steady state 
Nb = N/2. Atom- molecule oscillations are absent be- 
cause as soon as the atomic condensate is fully depleted, 
a(t) — 0, the nonlinear Rabi frequency in Eqs. [18] and 
[T51 vanishes leading to b — and a = 0. This behavior is 
the same as earlier mean field studies of atom-molecule 
conversion via Feshbach resonances or photo-association 
with undepleted lasers [H], H2, HH ■ In that case the an- 
alytic solution for the molecule number was shown to 
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be N b (t) = (N/2) tanh 2 (v / 2iV|fi|t) H2 where ft is atom- 
molecule coupling, which can be taken to be equivalent 
to a time independent g(e)/2 in our model. 

The AMF solution N b = N/2 is an unstable equilib- 
rium [3ll ] and fluctuations of the atom or molecular fields 
would excite the system out of this state leading to atom- 
molecule oscillations. The PAMMF equations incorpo- 
rate the vacuum fluctuations of the atomic-molecular 
fields in the term g(2N)^ 1 / 2 e in Eq. [15] Because of 
this term, even for 5 = Xa = Xb — when there is 
full conversion into molecules (L z (t) = +1), one still has 
L- = g(2N)~ 1 ^ 2 e. Figure 4 shows that when g(2N)~ 1 / 2 e 
is eliminated from Eq. [15] PAMMF and AMF equations 
produce virtually identical results [34j . 

Besides the molecule number, the AMF makes different 
predictions for the photon number than PAMMF and 
QSCP. For 5 = Xa = Xb = 0, the PAMMF and QSCP 
results indicate that the photon number approaches the 
steady state N e = while the AMF solution predicts that 
the photon number grows like N e — (2e/"f) 2 as shown 
in Fig. 5. This behavior is again attributable to the 
atom-molecule quantum fluctuations. When the AMF 
solution reaches the unstable equilibrium, N b — N/2, the 
equation for the cavity field reduces to that of an empty 
cavity, e = e — je/2. By contrast, the PAMMF mean 
field equations have the steady state solution L z = +1, 
e = 0, and L_ = (e/g)(2/N)~ 3 / 2 for the resonant case. 

For nonzero values of 5, x a , or Xb, the AMF equations 
do exhibit oscillations between atoms and molecules since 
the atom-molecule transition is detuned from perfect res- 
onance, which prevents full conversion into molecules 
from occurring. For example, for small \5\, PAMMF show 
higher frequency larger amplitude oscillations but as \5\ 
is increased, the agreement between AMF and PAMMF 
becomes increasingly better until the solutions are again 
indistinguishable for \5\ ^> gVN as seen in Fig. 3. 

These results indicate that for off-resonant atom- 
molecule conversion, all three methods agree. For res- 
onant transitions, the PAMMF and QSCP show atom- 
molecule Rabi oscillations due to the inclusion of the 
matter field vacuum fluctuations while the QSCP addi- 
tionally includes atom-molecule-photon correlations that 
modify the effective Rabi frequency. 



B. Weakly Driven Cavity, £<7 

For weak driving, e -C 7, reservoir noise and cavity 
vacuum fluctuations are both much stronger than the 
coherent driving. In this case, AMF and PAMMF equa- 
tions are expected to give an inaccurate description. In 
fact, for the limiting case e = 0, both mean field the- 
ories predict no dynamics at all since in the absence of 
photon fluctuations there is nothing to initiate the atom- 
molecule conversion. 

First we show a comparison of the PAMMF and QSCP 
equations for JV» 1 and e < 7 in Fig. 6. As one can 
see for the QSCP, the conversion of atoms into molecules 
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FIG. 3: (Color Online) (a) Molecule number and (b) intracav- 
ity photon numbers calculated using the PAMMF equations 
(solid lines) and AMF equations (dash dot lines) as a func- 
tion of time t (in units of 7 _1 ). In all cases N = 10 6 , g/j = 

4 x 10~ , and e/7 = 100. Here red lines are 8 = \ a = Xb — 0; 
blue lines are S/'y — 0.0001 and Xa = Xb — 0; green lines 
are 5/7 = 1 and Xa — Xb = 0- Note that in the final case 
(green lines), PAMMF and AMF are indistinguishable. Also, 
PAMMF for S/'y = 0.0001 (blue solid line) and PAMMF for 

5 = (red solid line) are indistinguishable. 
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FIG. 4: (Color Online) Molecule number calculated using the 
PAMMF equations without the term g(2N)~ 1/2 e in Eq. [15] 
(red solid line) and AMF equations (blue solid line). Here 
N = 10 6 , g/f = 10~ 5 , e/ 7 = 100, 5 = X a=Xb=0- The two 
solutions are indistinguishable. 



starts much earlier. This is easily understood if one solves 
the QSCP equations for L z (t) using perturbation theory. 
Solving for L z (t) for short times from the initial condition 
L z (0) = —1, one finds that 

L z {t) » -1 + (tg^N) 2 (l - 1/N)(1 - t</6), (26) 

which is correct to order t 3 even for arbitrary e. The 
cavity driving s only contributes to L z (t) at order t 4 in 
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FIG. 5: (Color Online) Steady state photon number as a 
function of e/7 calculated using the QSCP equations (solid 
lines), PAMMF equations (dotted lines), and AMF equations 
(dashed dot lines). Here N = 10 6 and 5/7 = 5 x 10" 5 . Blue 
lines are for 8 — Xa = Xb — and red lines are 5/7=1 and 

Xa = Xfc = 0. 



perturbation theory with the term 

L z (i)| e = ( 3 ViV) 2 (l-l/7V)| £ | 2 t 4 /4. (27) 

From which we see that the cavity fluctuations dominate 
the short time behavior and it is only for later times that 
the atoms feel the affect of the cavity driving. 



x 10 4 




1 2 3 4 5 6 7 



yt 

FIG. 6: (Color Online) (a) Molecule number and (b) intra- 
cavity photon numbers calculated using the PAMMF equa- 
tions (blue dash dot lines) and QSCP equations (red solid 
lines) as a function of time tf for the parameters N — 10 5 , 
g/j = 1.5 x 1(T 4 , 8 = X a = Xb = and e/7 = 0.01. 



In Figs. 7 and 8 we show solutions of the density ma- 
trix (Eq. |SJ), QSCP, and PAMMF solutions for initial 
atom numbers N a — 64 and e = and O.I7. One can 
see that for times long enough for nearly full conversion 
to occur, the QSCP and density matrix solutions show 
very good agreement while the PAMMF either predicts 



no dynamics (e = 0) or dynamics that start significantly 
later (e = O.I7). The molecule number, Nb, for both the 
PAMMF and QSCP solution approach a steady value 
that is below that of the density matrix, iV& = N/2. 
This is a consequence of the c-number approximation for 
L z , which is only valid for N 3> 1. Our simulations 
indicate that this error vanishes as N — > 00 and both 
PAMMF and QSCP indicate full molecule conversion on 
resonance. Unfortunately, our solutions of Eq. [5] are lim- 
ited to a maximum of 64 atoms and photons. 
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FIG. 7: (Color Online) (a) Molecule number and (b) intra- 
cavity photon numbers calculated using the density matrix 
(Eq. 5) (red solid lines), PAMMF equations (green dashed 
lines), and QSCP equations (blue dashed dot lines) as a func- 
tion of time <7 for the parameters iV = 64, 3/7 = 0.01, 
8 = Xa = Xb = and e/7 = 0. 



V. CONCLUSIONS 

Here we have compared several different approaches 
for analyzing the dynamics of intracavity photoassocia- 
tion. It has been shown that traditional mean field the- 
ory, which replaces annihilation/creation operators with 
c-numbers, fails to provide an accurate description of the 
dynamics since it ignores the quantum fluctuations of 
both the matter fields and the cavity field. We developed 
two new approaches that go beyond traditional mean 
field theory by introducing a different operator represen- 
tation for the Hamiltonian similar to angular momentum. 
When mean field theory is applied to the new operator 
representation, the resulting equations properly incorpo- 
rate vacuum fluctuations of the matter fields, which are 
necessary for atom-molecule Rabi oscillations to occur. 
The second approach, QSCP, which goes one step further 
by also incorporating quantum fluctuations and reservoir 
noise of the cavity field, provides an accurate description 
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FIG. 8: (Color Online) (a) Molecule number and (b) intra- 
cavity photon numbers calculated using the density matrix 
(Eq. 5) (red solid lines), PAMMF equations (green dashed 
lines), and QSCP equations (blue dashed dot lines) as a func- 
tion of time ty for the parameters N = 64, g/j — 0.01, 
S = Xa = Xb = and e/j = 0.1. 



of the short time behavior that agrees with exact nu- 
merical solutions of the density matrix equations. This 
final method that incorporates photon noise works even 
when the cavity in not driven, which is when both mean 
field theories completely fail. Even more important, the 
QSCP equations can be easily solved numerically for ar- 
bitrary numbers of atoms and photons while the density 
matrix can only be solved for relatively small numbers of 
atoms and photons, typically less than a 100 each. 



In a future work we plan to improve upon the QSCP 
method in a manner that will allow us calculate both 
squeezing and number fluctuations of the fields. 



The work was supported by the National Science Foun- 
dation award no. 0757933. 
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